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Abstract 

In this paper, a proposal for the restriction of the Euclidean functional integral to a region 
free of infinitesimal Gribov copies in linear covariant gauges is discussed. An effective action, 
akin to the Gribov-Zwanziger action of the Landau gauge, is obtained which implements the 
aforementioned restriction. Although originally non-local, this action can be cast in local 
form by introducing auxiliary fields. As in the case of the Landau gauge, dimension two 
condensates are generated at the quantum level, giving rise to a refinement of the action 
which is employed to obtain the tree-level gluon propagator in linear covariant gauges. A 
comparison of our results with those available from numerical lattice simulations is also 
provided. 


1 Introduction 

Analytical approaches to Quantum Chromodynamics (QCD) face many problems due to the 
fact that at low energy scales, where many important physical phenomena as confinement and 
chiral symmetry breaking take place, perturbation theory breaks down. Moreover, a complete 
non-perturbative framework is not yet at our disposal. On the other hand, lattice numerical 
methods have provided firm evidence for those non-perturbative effects. Let us underline that 
the interplay between numerical and analytical analysis has been proven to be very fruitful over 
the last decade. 

An analytical approach which is receiving growing attention is the Gribov-Zwanziger set up. 
In this framework, a suitable modification of the Faddeev-Popov path integral quantization 
procedure is worked out and non-trivial non-perturbative results emerge. As first noticed by 
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Gribov [lj, after gauge fixing, one still has gauge field configurations which satisfy the gauge 
condition and are connected to each other through a gauge transformation: the so-called Gribov 
copies. This implies that a residual redundancy is still present in the Faddeev-Popov quantization 
formula. Although this fact was originally pointed out in Landau and Coulomb gauges |1], it was 
proven by Singer [2] that this is not a peculiarity of a particular gauge choice, being a general 
feature of the gauge fixing procedure. 

In [Tlpl a partial solution to this problem was proposed in the Landau gauge, = 0, for which 

the corresponding Faddeev-Popov operator M ab , 

M ab = -6 ab &+gf abc A°d li , with = 0, (1) 

is hermitian M, so that its eigenvalues are real. The original proposal by Gribov |I] in order to 
take into account the existence of gauge copies was that of restricting the domain of integration 
in the functional Euclidean integral to a certain region Q in field space defined as the set of 
all field configurations fulfilling the Landau gauge and for which the Faddeev-Popov operator is 
strictly positive, namely 

= 0, M ab > 0 } . (2) 

The region Q is called the Gribov region. It enjoys several properties [5j: 1) it is a bounded 
region in field space, 2) it is convex, 3) all gauge orbits cross Q at least once. In particular, this 
last property ensures that a gauge field configuration located outside of the Gribov region is a 
copy of a configuration located inside f 1. This important result gives a well defined supportQ to 
Gribov’s original proposal of restricting the functional integral to the region fi. 

Therefore, according to Gribov’s proposal p], for the partition function of quantized Yang-Mills 
action in the Landau gauge we write 

Z = [ [T>$] e -( 5YM+S 8f) , (3) 

Jn 

with (Sym + S'gf) given by eq. (l65l) with a = 0. Later on, Zwanziger [6] was able to show that the 
restriction of the domain of integration to the region tt can be effectively implemented by adding 
to the starting action an additional non-local term Hl(A), known as the horizon function. More 


precisely 

/ [D$] e - ( SYM+s e f ) = / [D$]e _s G Z) 

Jn J 

(4) 

where 

Sh z = Sym + Sgf + 1 A H l (A) - 4 V^(N 2 - 1), 

(5) 

and 

H L {A) =g 2 f d A x f abc A\\ f dec Al , 

(6) 

*See [3ia for a pedagogical introduction to the Gribov problem. 

^Nowadays, it is known that the Gribov region P is not completely free from Gribov copies, see for example 
[4j and references therein. There are still additional copies inside P. A smaller region, known as the fundamental 
modular region (FMR) has been identified which is completely free from Gribov copies. Though, our present 
knowledge of the FMR is not as good as our knowledge of the Gribov region P, which will be the main subject 


of the present investigation. 
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with [Ad -1 ] denoting the inverse of the Faddeev-Popov operator, eq.(0). The mass parameter 7 
appearing in expression ([5]) is known as the Gribov parameter. It is not a free parameter, being 
determined in a self-consistent way by the gap equation | 6 J 

(H l ) = 4V(N 2 - 1). (7) 

where the vacuum expectation values (Hl) has to be evaluated with the measure defined in 
eq.®). 

The action in eq.([5]) is called the Gribov-Zwanziger (GZ) action. It is a non-local action due to 
the presence of the horizon function. However, it can be cast in local form [4j by the introduction 
of a pair of commuting fields ((( b ) and a pair of anti-commuting ones (<7“ 6 , namely 


Sh z = Sym + ,%- 

+ 7 2 f d A x gf abc A“(<p + - 4 7 4 F(IV 2 - 1). (8) 

Since the auxiliary fields (kp ab , and (uj a ^,uj ab ) appear quadratically in expression ((SI), they 
can be easily integrated out, giving back the non-local expression Q. Remarkably, the local 
action (J5j) turns out to be renormalizable to all orders HIE!, while implementing the restriction 
to the Gribov region fl in the Landau gauge. 

The Gribov-Zwanziger action in the Landau gauge is plagued by dynamical non-perturbative 
instabilities and the generation of dimension two condensates is favoured 0019]. These con¬ 
densates can be taken into account by further modifying the original action, giving rise to the 
so-called refined Gribov-Zwanziger (RGZ) action (7. 8][9 . given by 

Ska z = Ski + ^ / d‘x KK ~ V 2 / d'x (« - «) . P) 

As much as the Gribov parameter 7 2 , the parameters (m,M) are not free, being generated 
dynamically in a self-consistent way, see EM- 

The RGZ action describes gluon and ghost propagators in very good agreement with the most 
recent lattice results [ 1(3 mm- In particular, for the the tree level gluon propagator one 
obtains 


(A«(k)Al(-k)) L 



k tl kA k 2 + M 2 

k 2 ) ( k 2 + m 2 )(k 2 + M 2 ) + 2Ng 2 7 4 ' 


( 10 ) 


Also, the RGZ framework gives results for glueball masses which compare well with numerical 
results, see mm and references therein. 

It is worth underlining here that the construction of the GZ and RGZ actions in the Landau 
gauge follows directly from the excellent understanding of the properties of the Gribov region kl 
of the Landau gauge which has been achieved so far, as mentioned before. Willing to address the 
Gribov issue in other gauges, the first and necessary step towards the construction of a suitable 
non-perturbative action is the study of the corresponding Gribov region, whose properties are 
encoded in the structure of the Faddeev-Popov operator associated to the chosen gauge. This 
is, for example, the way in which the Gribov issue has been addressed in the maximal Abelian 
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gauge [ei ces cm nsi ns], where it has been possible to characterize a few properties of the 
corresponding Gribov region, paving the construction of the analogue of the RGZ action in this 
gauge. Let us mention that, besides the Landau and maximal Abelian gauges, a successful 
characterization of the Gribov region and of the ensuing GZ and RGZ actions has been worked 
out also in the non-covariant Coulomb gauge [20] EU [ 22 , 22 j [23, 24]. All these three gauges 
share the very important feature that the corresponding Faddeev-Popov operators are hermitian. 
Thus, they have real eingevalues, a property which makes the definition of the corresponding 
Gribov regions very clear and transparent. In addition, all these gauges can be defined through 
the minimization of suitable auxiliary functionals. For example, the Landau gauge, = 0, 

can be obtained as the stationary condition with respect to the gauge transformations of the 
auxiliary functional f d 4 xA^A^. Moreover, the second variation of f d 4 xA“A“ with respect to 
the gauge transformations yields precisely the Faddeev-Popov operator of the Landau gauge, 
see (Tj. As a consequence, the Gribov region D of the Landau gauge, eq.©, can be defined as 
the set of all relative minima of the auxiliary functional f d 4 xA^A^ in field space. A similar 
feature holds in both maximal Abelian and Coulomb gauges. This important property has 
deep mathematical consequences and allows for a lattice formulation of these gauges, whose 
predictions can be thus tested from the numerical side. 

In general, addressing the issue of the Gribov copies in an arbitrary gauge is a very difficult task, 
as many properties of the Landau, Coulomb and maximal Abelian gauges are lost. This is, for 
instance, the case of the linear covariant gauges which are the object of the present investigation. 
Here, unlike the Landau, Coulomb and maximal Abelian gauges, the Faddeev-Popov operator 
lacks hermiticity and, moreover, no auxiliary functional is at our disposal, making the treatment 
of the Gribov copies highly non-trivial. Despite these difficulties, a first attempt to address the 
Gribov problem in the linear covariant gauges was discussed in [25] , under the assumptions that 
the gauge parameter a present in these gauges was considered to be infinitesimal, i.e. a < Cl. 
This hypothesis allowed the authors [25] to introduce a region in field space which is free from 
infinitesimal Gribov copies. More precisely, according to [25], this region is defined by requiring 
that the transverse component of the gauge fielcJH = (S^—d fl d l ,/d 2 )A^ belongs to the Gribov 

region H of the Landau gauge, eq. (f 2 l) . and that the longitudinal component = ( d^d u /d 2 )A “ 
remains free. 

In this work we pursue the analysis of the linear covariant gauges started in [25]. In particular, we 
identify a suitable region Hlcg hi field space which turns out to be free from infinitesimal copies, 
without relying on the assumption a •Cl made in [25] • More precisely, as we shall see in the next 
section, we will be able to eliminate all infinitesimal copies related to zero modes of the Faddeev- 
Popov operator which are smooth functions of the gauge parameter a. Also, the restriction of the 
domain of integration in the Euclidean functional integral is effectively implemented through the 
introduction of the correspondent non-local horizon function which can be localized by means 
of a set of auxiliary fields. Though, the localization procedure in the linear covariant gauges 
displays new features with respect to the Landau gauge. Therefore, a local Gribov-Zwanziger 
action which eliminates infinitesimal Gribov copies in linear covariant gauges is established. The 
gap equation which determines the Gribov parameter 7 is analyzed at one-loop order, with the 
important physical outcome that no gauge parameter dependence emerges for 7 . As in the case 
of the Landau, maximal Abelian and Coulomb gauges, this action is plagued by further non- 
perturbative dynamical effects related to dimension two condensates, resulting in a refinement 

^We remind that the class of the linear covariant gauges is defined by the condition: — ab a = 0, where 

a is a gauge parameter. Unlike the Landau gauge, the gauge field A“ has now both transverse and longitudinal 
components, i.e. A“ = A“ T + 
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of the starting Gribov-Zwanziger action. The tree-level gluon propagator stemming from the 
refined Gribov-Zwanziger action in linear covariant gauges is evaluated and compared with the 
most recent numerical lattice lattice data. 

The work is organized as follows. In Sect J2] we elaborate on the definition of the region Qlcg 
which allows to eliminate the infinitesimal copies in the linear covariant gauges. In Sectj3]we 
implement the restriction of the domain of integration in the functional integral to the region 
Ulcg by introducing a suitable non-local horizon function. In Sect|4]we show how to localize 
the action defined in the previous section by introducing a set of auxiliary fields. The resulting 
action is identified with the local Gribov-Zwanziger action in linear covariant gauges. In SectJSl 
the BRST soft breaking of the Gribov-Zwanziger action is briefly analyzed. A consistency 
check is realized in SectlHJ where we show that the gap equation which determines the Gribov 
parameter 7 turns out to be independent from the gauge parameter at one-loop order. Further, 
the origin of dimension two condensates by dynamical effects is analyzed in Sect|7]as well as the 
refinement of the Gribov-Zwanziger action. In SectJH] we present the tree-level gluon propagator 
and compare it with the available results from numerical lattice simulations. Finally, we present 
our conclusions and perspectives. 


2 Proposal for a Gribov region in linear covariant gauges 

As discussed in Sect 03 a partial resolution of the Gribov problem in the Landau gauge requires 
the restriction of the domain of integration in the path integral to the so-called Gribov region 
P, eq. (J23). In the Landau gauge, since the Faddeev-Popov operator is hermitian, the definition 
of 0 as the region in field space where the Faddeev-Popov operator is positive has a clear and 
well defined geometrical meaning, as already mentioned before. This is not the case of linear 
covariant gauges, where the gauge condition reads 

- ab a = 0 , ( 11 ) 

where a. is an arbitrary gauge parameter and b a a Lagrange multiplier held. This gauge condition 
coincides with Landau condition when a = 0. However, for values of a different from zero, the 
gauge held acquires a longitudinal component, which is absent in the Landau gauge. If we 
perform an inhnitesimal gauge transformation on ea. llll) and impose the gauge condition (HID 
on the gauge transformed held, we obtain 


M°\ 6 = -d„Dfi b = 0, (12) 

M ab = -d.Df = + gf acb AD = -S ab d 2 - agf ach b c - gf^A^ , (13) 

with £ a being the inhnitesimal parameter of the gauge transformation. The operator JA ab in 
eq. ITdl) is the Faddeev-Popov operator of the linear covariant gauges and eq. (fT2l) is the corre¬ 
sponding zero mode equation giving rise to inhnitesimal Gribov copies, namely 

A; = A; + D ; b e, d„A“=d,A*=ab a . (M) 

Unlike the Faddeev-Popov operator of the Landau gauge, the operator M ab in eq. m is not 
hermitian. This fact makes the characterization of the analogue of the Gribov region in the 
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linear covariant gauges highly non-trivial. We highlight that the elimination of the Gribov 
copies using the approach followed by Gribov [lj and by Zwanziger Eli requires the use of a 
hermitian operator. The understanding of the properties of the so-called Gribov region is crucial 
in order to deal with the problem of the gauge copies. As already mentioned, for a generic choice 
of the gauge condition, the properties of the corresponding Gribov region are encoded in the 
associated Faddeev-Popov operator and no systematic way of defining it is known. 

Nevertheless, an attempt to define a region free from infinitesimal copies in linear covariant 
gauges was done in [25] . In this work, a region free from infinitesimal copies was identified, 
under the assumption that the gauge parameter a was infinitesimal, i.e. a < 1, The region 
introduced in [25j is defined by demanding that the transverse component of the gauge field, 
A“, belongs to the Gribov region 12 of Landau gauge, eq. ([2]). Here, we 
extend this result for a finite value of a. The extension relies on the following theorem. 

Theorem 2.1. If the transverse component of the gauge field A^ a = (V - € Q, then 

the equation = 0 has only the trivial solution £ b = 0. 

Proof. Since A°^ £ H by assumption, the operator 

M Tab _ _ § abQ 2 + gf abc A TcQ^ , (i5) 

is positive definite and, therefore, is invertible. As a consequence, ea. (H2l) can be rewritten as 


C(x,a) = -g[{U T )~Y 

= -ga [(M T ) -1 ] f dbc 8 „ ( b c ) e 6 ) , (16) 

where the gauge condition CD was used. We consider here zero modes f(x,a) which are 
smooth functions of the gauge parameter a. This requirement is motivated by the physical 
consideration that the quantity f(x, a ) corresponds to the parameter of an infinitesimal gauge 
transformation. On physical grounds, we expect thus a regular behaviour of £(x,a) as function 
of a , i.e. infinitesimal modifications on the value of a should not produce a drastic singular 
behaviour of £, a feature also supported by the important fact that an acceptable zero mode has 
to be a square-integrable function, i.e. f d 4 x f a f a < oo. Also, the a-dependence should be such 
that in the limit a —> 0 we recover the zero-modes of the Landau gauge. Therefore, we require 
smoothness of £(x, a) with respect to a. Thus, for a certain radius of convergence R, we can 
write the zero-mode £(x, a) as a Taylor expansion in a, 


$ a (x,a) = 52 a n {Z(x). 

71=0 

For such radius of convergence, we can plug ea. (fT71) into ea. (fT6l) . which gives 

OO OO / / rj \ 

= -X > n+1 [(m T )~ Yf dbc d , ((^ c ) & 

n =0 n =0 ^ ^ 



(17) 


(18) 
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Since a is arbitrary, eg. (1181) should hold order by order in a, which implies 

Co =0 => 4>o = -a [Pvt T )- 1 ]“'/ <fc c>„ ((I? 60 ) «o) = 0. (19) 

at zeroth order. Therefore, 

a = O0S = 0 , (20) 

and by recursion, 

c = c<_i = 0, Vn. (21) 

Hence, the zero-mode £(x,a) must be identically zero within R. Due to the requirement of 
smoothness, i.e. of differentiability and continuity of £, the zero-mode must vanish everywhere. 

□ 

A comment is in order here. As underlined before, Theorem m holds for zero-modes which 
have a Taylor expansion in powers of a. Although this does not seem to be a very strong 
requirement, since we expect smooth functions of ct, i.e. small perturbations on a should not 
result on abrupt changes on £, one could think about the possibility to have zero-modes which 
might eventually display a pathological behaviour, i.e. which could be singular for some values 
of the gauge parameter a. For this reason, we shall refer to the zero-modes that admit a Taylor 
expansion as regular zero-modes. 

Motivated by the previous Theorem, we introduce now the following Gribov region Dlcg i n the 
linear covariant gauges: 

Definition 2.1. The Gribov region Olcg in linear covariant gauges is given by 

^lcg = {tL“, - ab a = 0 , M Tab > o} , ( 22 ) 

where the operator M. Tab is given in eq. (H5l) 

From the previous Theorem, it follows that the region Hlcg is f ree from infinitesimal Gribov 
copies which are regular, i.e. smooth functions of the gauge parameter a. 

Although the understanding of the Gribov issue in the linear covariant gauges cannot certainly 
be compared to that achieved in the Landau gauge, the possibility of introducing the region 
Hlcg which is free from infinitesimal copies can be regarded a first important step in order 
to face this highly non-trivial problem. Let us now focus on the construction of the horizon 
function in this class of gauges 


3 The horizon function and the Gribov-Zwanziger action in lin¬ 
ear covariant gauges 

Definition 12.11 provides a consistent candidate for the Gribov region in linear covariant gauges. 
Following the approach employed by Gribov in [Tj and generalized by Zwanziger [fij, we should 
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restrict the path integral to the region ^lcG) Le. 


Z= [£$] e" (SYM+5 sf) , (23) 

-'^LCG 

where $ represents all fields of the theory, see Appendix. [A] for the conventions used. From 
equation ([22]) . one immediately sees that the region flLCG is defined by the positivity of the 
operator M T which contains only the transverse component of the gauge field, ea. (115l) . In other 
words, M t is nothing but the Faddeev-Popov operator of the Landau gauge. As a consequence, 
the whole procedure performed by Gribov [1] and Zwanziger [5] in the case of the Landau gauge, 
can be repeated here, although one has to keep in mind that the restriction of the domain of 
integration to the region Ulcg affects only the transverse component of the gauge field, while 
the longitudinal sector remains unmodified. Therefore, following EE], for the restriction to the 
region Ulcg we write 


/ [D$] Q~( s ™+S S i) = / p$] e -S GZ ; (24) 

•I^LCG J 

where the action Sq z is given by 

5gz = Sym + % + 7 4 iL(A T ) - - 1). (25) 

The quantity H(A T ) is the non-local horizon function which depends only on the transverse 
component of the gauge field A T , namely 

H(A t ) =g 2 J d 4 x / adc A^ c [(M T ) _1 ] a6 / Me ^ e ■ (26) 

The parameter 7 in ea. (125l) is the Gribov parameter. As in the case of the Landau gauge, it is 
determined in a self-consistent way by the gap equation 

(. H(A t )) = 4V(N 2 - 1), (27) 

where the vacuum expectation value ( H(A T )) has to be evaluated now with the measure defined 
in eg. f[24l) . 

The effective action (1251) implements the restriction to the region Ulcg- Here, an important 
feature has to be pointed out. Formally, the horizon function (|26l) is the same as in Landau 
gauge. However, in this case, although the longitudinal component of the gauge field does not 
enter the horizon function, it appears explicitly in the action Sgz- As we shall see, this property 
will give rise to several differences with respect to the Landau gauge. Another important point 
to be underlined concerns the vacuum term 4V"y 4 (N 2 — 1) in expression (125[) . This term is 
related to the spectrum of the operator JA T [3] which does not depend on the gauge parameter 
a. Therefore, at least at the level of the construction of the effective action S'gz, eq. ([25]h the 
vacuum term is independent of a. Of course, we should check out if quantum corrections might 
eventually introduce some dependence from a in the vacuum term. This would require a lengthy 
analysis of the renormalizability properties of 5 gz ; which is under current investigation [26] ■ 

From now on, we will refer to the action (1251) as the Gribov-Zwanziger action in linear covariant 
gauges. As it happens in the Landau gauge, this action is non-local, due to the non-locality of 


the horizon function. However, as we shall see in the next section, it is possible to localize this 
action by the introduction of a suitable set of auxiliary fields. Here, differences with respect to 
the case of the Landau gauge will show up, due to the unavoidable presence of the longitudinal 
component of the gauge field. 

4 Localization of the Gribov-Zwanziger action in linear covari¬ 
ant gauges 

In order to have a suitable framework to apply the usual tools of quantum field theory, we have 
to express the action (1251) in local form. In the case of linear covariant gauges, the localization is 
not as direct as in Landau gauge. The difficulty relies on the fact that the horizon function (|26l) 
has two kinds of non-localities. The first one is the same as in Landau gauge, i.e. the presence 
of the inverse of the Faddeev-Popov operator . The other one follows from the fact that the 
decomposition of the gauge field into transverse and longitudinal components is non-local, see 
eq. dSzD. Therefore, if we apply the same procedure used in the Landau gauge, the localization 
of the horizon function would give rise to a term of the type 

J d A x 91 2 r bc A T ^ + (p)f, (28) 

which is still a non-local term, due to the presence of the transverse component A“ T . However, 
it is possible to localize the action ([25]) using an additional step. First, let us write the transverse 
component of the gauge field as 



K~h 


a 

M ’ 


(29) 


where the field /i“ will be identified with the longitudinal component, i.e. we shall impose that 


h 


a 


f) f) 

u H u v a a 

d 2 


(30) 


In other words, we introduce an extra field and state that the transverse part of the gauge 
field can be written in a local way using ea. (l29D . Clearly, we must impose a constraint to ensure 
that, on-shell, eq. (l29l) is equivalent to the usual decomposition (fBD) . Before introducing this 
constraint, we rewrite the horizon function in terms of As a matter of notation, we will 
denote M r as 3Vt(A — h), when the transverse gauge field is expressed in terms of /i“. The 
horizon function (1261) is now written as 


H(A, h) = g 2 j d 4 x f adc (A ]= - ^)[(M(A - f Me ^ - h %). (31) 

The constraint given by ea. (l30l) is imposed by means of a Lagrange multiplier A“, he., by the 
introduction of the term 


Sx 


/ 


d 4 x A a Jd 2 K - d lt d v A a v ) 


(32) 
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Also, a second constraint is required in order to ensure that the held (A“ — h“) is transverse, 


8Ml~K) = o, 

which is implemented by the following term 


( 33 ) 


S v = J d*xv a d^A“-h“), (34) 

with v a being a Lagrange multiplier. Therefore, the introduction of the extra held /i“ in eq. (1.31 h 
by means of the constraints (1301) and (1331) , implemented by the terms (1321) and (fMl) , provides an 
action S GZ 


S' GZ = S YM + % + S X + S V + 7 4 H(A, h) - 4R 7 4 (1V 2 - 1) (35) 

which is on-shell equivalent to the the non-local Gribov-Zwanziger action (1251) . The introduction 
of the helds /i“, A“ and v a has to be done through a BRSlI doublet [27] to avoid the appearance 
of such helds in the non-trivial part of the cohomology of the BRST operator s , a property which 
will be important for the renormalizability analysis, [26). Therefore, we introduce the BRST 
doublets (/*“,£“), (A“,A“) and (r a ,u a ), i.e. 

= S A“ = 0 

= 0 , ST a = V a 

sA“ = A“, su a = 0 (36) 

and dehne the following BRST exact terms 

S \\ = s J d 4 x A “(d 2 h; - d,d v Al) = J d A x A ;(d 2 h“ - d^Al) , (37) 

- J dW^dX + d^D^) 

Stv = s J d A x r a (9 /i (A“ - h^) = J d 4 x 'i/’d^A" - /i“) 

+ j d A XT a d^Dfc b + £“). (38) 


The terms (j37[) and (1381) implement the constraints (1301) and ()33[) in a manifest BRST invariant 
way. What remains now is the localization of the term (1311) . Since this term has just the usual 
non-locality of the Gribov-Zwanziger action in the Landau gauge, given by the inverse of the 
operator M(A — h), we can localize it by the introduction of the same set of auxiliary helds 
employed in the case of the localization of the horizon function in the Landau gauge [1]. Thus, 
the term (1311) is replaced by the the following local expression 

§ For the usual BRST invariance of the Faddeev-Popov action in linear covariant gauges we remind the reader 
to Appendix m 
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S H = s J d 4 x w“ c M ab (A - h )^ c + 7 2 5 J d 4 x f abc (A“ - h\\){<p + ^ 

= J d 4 x (q%d v D<?<p* - uY?d v D?u$ - gf adb {d^){Dt e c e )^) 

+ J d 4 x - gf adb {d^)ht^ - gf adb {d^)ii^) • 

( 39 ) 

where 


«*>?=<, 
»< = 0 , 


» 7 * = ft 
s?t = 0 . 


Finally, the action Sgz given by 


(40) 


Sgz = Sym + -5gf + Si x + S TV + Sh (41) 

is local and, on-shell, equivalent to the non-local action ([26]) . Explicitly, Sqz is written as 


Sgz = \ j d 4 x + j d 4 x b a (d^A* - % a ) + j d 4 x c a d,D ab c b 

+ J d 4 x A“(<9 2 /i“ - dfAAl) + J d 4 x A“(<9 2 £“ + d^d u D“ h c b ) 

+ J d 4 x v a d^Al - h%) + J d 4 x r a d^{D ab c b + £“) + J d 4 x (v^d v D^<p^ 

- - gf^id^Xl^tf)^) + j d 4 x ( gf adb {d v ^)h d u 

- gr db {d u u a ;)h d u ^ - S r(^<)6f) + 7 2 9 J d 4 x f abc (A “ - /£)(*, + $>)£ 

- 4E 7 4 (A ^ 2 - 1) , (42) 

and we will refer to it as the local Gribov-Zwanziger action in linear covariant gauges. We 
underline that, in the limit 7 —>■ 0, the term (|39D can be trivially integrated out to give a 
unity. The remaining action is simply the gauge fixed Yang-Mills action with the addition of 
the constraints over h“. These constraints are also easily integrated out, so that the resulting 
action is simply the usual Yang-Mills action in linear covariant gauges, see Appendix [Aj 

Let us end this section by noticing that, in the local formulation, the gap equation (1271) takes 
the following expression 


dS 1 = Q 

d 7 2 


(43) 


11 


where £ 7 denotes the vacuum energy, obtained from 


e~ v£l = J [D$] e~ Scz . (44) 

with 4> being the complete set of fields, he., the usual ones from the Faddeev-Popov quantization 
and the auxiliary fields introduced to implement the constraints and to localize the Gribov- 
Zwanziger action. 

5 BRST soft breaking of the Gribov-Zwanziger action in linear 
covariant gauges 


As it happens in the case of the Gribov-Zwanziger action in the Landau and maximal Abelian 
gauges eq 0 la m, ii3 im uni m ns nn m m, gsi m mg, the expression (1421) is not invariant 
under the BRST transformations () 66 l) . (1361) and (1401) . The only term of the action which is not 
invariant under BRST transformations is 

57 2 j d 4 x f abc (A“ - /*“)(</> + , (45) 


giving 


•sSgz = A 7 2 = gY 


1 


d A x f 


abc 


~(D?c d + £)(*> + vt + K * KX 


(46) 


From eq. (!46p . we see that the BRST breaking is soft, i.e. it is of dimension two in the quantum 
fields. This is precisely the same situation of the Landau and maximal Abelian gauges. The 
restriction of the domain of integration in the path integral to the Gribov region generates a 
soft breaking of the BRST symmetry which turns out to be proportional to the parameter y 2 . 
As discussed before, when we take the limit 7 —> 0, we obtain the usual Faddeev-Popov gauge 
fixed Yang-Mills action which is BRST invariant. Thus, the breaking of the BRST symmetry 
is a direct consequence of the restriction of the path integral to the Gribov region. We also 
emphasize that, although the gauge condition we are dealing with contains a gauge parameter 
a, the BRST breaking term does not depend on such parameter, due to the fact that the horizon 
function takes into account only the transverse component of the gauge fields, as eq. (1261) shows. 


6 Analysis of the gap equation at one-loop order 

As discussed above, in the construction of the effective action which takes into account the pres¬ 
ence of infinitesimal Gribov copies a non-perturbative parameter 7 , i.e. the Gribov parameter, 
shows up in the theory. However, this parameter is not free, being determined by the gap equa¬ 
tion (1431) . In the Landau gauge, the Gribov parameter encodes the restriction of the domain of 
integration to the Gribov region P. Also, physical quantities like the glueball masses were com¬ 
puted in the Landau gauge na nn, exhibiting an explicit dependence from 7 . It is therefore of 
primary importance to look at the Gribov parameter in linear covariant gauges, where both the 
longitudinal component of the gauge field and the gauge parameter a are present in the explicit 
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loop computations. We should check out the possible (in)dependence of 7 from a. Intuitively, 
from our construction, we would expect that the Gribov parameter would be independent from 
a, as a consequence of the fact that we are imposing a restriction of the domain of integration 
in the path integral which affects essentially only the transverse sector of the theory. Moreover, 
the independence from a of the Gribov parameter would also imply that physical quantities like 
the glueball masses would, as expected, be ct-independent. To obtain some computational con¬ 
firmation of the possible a-independence of the Gribov parameter, we provide here the explicit 
computation of the gap equation at one-loop order. 

According to eqs. f[43]) , ([141) . the one-loop vacuum energy can be computed by retaining the 
quadratic part of Sqz and integrating over the auxiliary fields, being given by 

—yA 1 ) 
e 7 = 


[' DA] e ^ 


(p) AM(-p)+4 C 7 4 (AT a --1) 


(47) 


where 


7\ab _ ra6 




2 2Ng 2 ^ 

p + 


+ P^Pv 


1 — a 


2JV 5 2 7 4 


P J \ V a / P 

Performing the functional integral over the gauge fields, we obtain 


(48) 


V£^ = ^Tr InAj - 4U 7 4 (IV 2 - 1). 


(49) 


The remaining step is to compute the functional trace in ea. (l49IK This is a standard computation, 
see @]. Using dimensional regularization, we have 


,( 1 ) {N 2 — l){d — 1) 

7 2 




2Ng 2 ^ \ 

p 2 J 


d^(N 2 - 1), 


(50) 


where d = 4 — e. We see thus from ea. (|50l) that the one-loop vacuum energy does not depend 
on a and the gap equation which determines the Gribov parameter is written as 


d£? (' d-l)Ng 2 f d d p 1 

d'y 2 d J (2^)^ p 4 + 2A^5 2 7 4 

This equation states that, at one-loop order, the Gribov parameter 7 is independent of a and, 
therefore, is the same as in the Landau gauge, which agrees with our expectation. Although 
being a useful check of our framework, it is important to state that this result has to be extended 
at higher orders, a non-trivial topic which is already under investigation [26]. 


7 Dynamical generation of dimension two condensates 

The inclusion of dimension two condensates, see M, is an important aspect of the Gribov- 
Zwanziger set up in the Landau gaug^H, giving rise to the so-called RGZ action whose predictions 

"it is worth underlining that these condensates are also present in the maximal Abelian and Coulomb gauges, 

see ng esi ezi iia usual. 


13 















are in very good agreement with the most recent lattice numerical data. We underline that the 
inclusion of such condensates is not done by hand. They emerge as non-perrurbative dynamical 
effects due to a quantum instability of the Gribov-Zwanziger theory. In this sense, the RGZ 
action M takes into account the existence of these condensates in an effective way already 
at the level of the starting action. Here, we expect that, in analogy with the Landau gauge, 
dimension two condensates will show up in a similar way. In fact, the presence of these dimen¬ 
sion two condensates can be established in a very simple way through a one-loop elementary 
computation, which shows that the following dimension two condensates 

, («-«■>. (52) 

are non-vanishing already at one-loop order, being proportional to the Gribov parameter 7 . In 
particular, it should be observed that the condensate (A^ a A'^ a ) contains only the transverse 
component of the gauge field. This is a direct consequence of the fact that the horizon function 
of the linear covariant gauges, eq. (l26l) . depends only on the transverse component A °^. In order 
to evaluate the condensates ( A T A T ) and (<pip — Colo) at one-loop order, one needs the quadratic 
part of the Gribov-Zwanziger action in linear covariant gauges, namely 


o( 2 ) _ 

D GZ — 


/ 


d 4 x 


a 


-Al5, L J ab d 2 Al + ( 1 -) A^d.AAl) + ?5 ab d 2 c b 




1 -ac ra6 o2, ^6c -acea6o2, ,6c , '2 „ / j 4 „ rabc a a I c 

+ 0 O ip-u„d d uo + 7 9 / a x / A v I o^ u — 






/■ 


pabc a a 


d^d u 

d 2 


(<P + <P) 


-z\bc 


- 4H 7 4 (A^ 2 - 1). 


(53) 


Further, we introduce the operators f d 4 xA T A T and f d 4 x(cpp — coco) in the action by coupling 
them to two constant sources J and m, and we define the vacuum functional £(m, J) defined by 

e -V£(m,J) = 


[Dd>] q~~ s gz +j f d4x (K c K c -^ Cu Z c )- m f d4x A z{ 5 ^- C - 


(54) 


It is apparent to check that the condensates (A T A T ) and {tptp — oooo) are obtained by differenti¬ 
ating £(m, J ) with respect to the sources (J, m), which are set to zero at the end, i.e. 


/ -ac.jae -ac, ,ac\ 

VPu. Vu -U.LO } = - 


d£(m , J) 


(K a Al a ) 


dJ 
d£(m , J) 


J=m =0 


dm 


J=m =0 


A direct computation shows that 


e -V£(m,J) = e -^Tr lnA^+4y 7 4 (7V 2 -l) 


(55) 


(56) 


with 


A ab 

flV 


5 ab 




2 -f 4 g 2 N 
k 2 + J 


+ 2 m 


+ k^kp 


1 — a 
a 


2'y 4 g 2 N 2 m\ 

k 2 (k 2 + J) ~T 2 ) 


(57) 
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Evaluating the trace, we obtaiiJI 


T . (d — 1)(N 2 — 1) f d d k , / )2 2 7 4 5 2 1V , , 4/at2 

£(m, J ) =--- / , In ( fc 2 + —-— + 2m ) - d^(N 2 - 1) 


(2tt) 


k 2 + J 


(58) 


Eg. (1551) and (1551) gives thus 

- «> = 7 VN(N 2 - 1 )(d -1) j 


d d k 1 


(27r) d k 2 (fc 4 + 2g 2r j 4 N) 


(59) 


and 


(AjMj a ) = - 7 4 (JV 2 -l)(d-l) 


d d k 1 2 g 2 ./V 

(27r) d k 2 (fc 4 + 2g 2 'y 4 N) ’ 


(60) 


where we have employed dimensional regularization. Eq. (|59l) and eq. dBOl) show that, already at 
one-loop order, both condensates (A T A T ) and (<^</? — uu) are non-vanishing and proportional to 
the Gribov parameter 7 . Notice also that both integrals in eqs. (|59l) . (160l) are perfectly convergent 
in the ultraviolet region by power counting. We see thus that, in perfect analogy with the case of 
the Landau gauge, dimension two condensates are automatically generated by the restriction of 
the domain of integration to the Gribov region, as encoded in the parameter 7 . As shown in Ell, 
the presence of these condensates can be taken into account directly in the starting action giving 
rise to the refinement of the Gribov-Zwanziegr action. Also, higher order contributions can be 
systematically evaluated through the calculation of the effective potential for the corresponding 
dimension two operators by means of the Local Composite Operator technique, see 0136]. 

In the present case, for the refined version of the Gribov-Zwanziger action which takes into 
account the presence of the dimension two condensates, we get 


Srgz = Sgz + ^J d 4 x (a; - h“)(A“ - hi) - M 2 J d±x (<^V“ b - «) , (61) 

where the parameters (m, M) can be determined order by order in a self-content way through the 
evaluation of the corresponding effective potential, as outlined in the case of the Landau gauge 
[9]. Let us remark here that the calculation of the vacuum functional £(m, J ), eq. (15^1) . done in 
the previous section shows that, at one-loop order, these parameters turn out to be independent 
from a. The study of the effective potential for the dimension two operators f d 4 xA T A T and 
f d 4 x(ifip — Cdu) will be thus of utmost importance in order to extend this feature to higher 
orders. 

We are now ready to evaluate the tree level gluon propagator in the linear covariant gauges. 
This will be the topic of the next section. 

7n this computation we are concerned just with the contribution associated to the restriction of the path 
integral to Ulcg ■ 
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8 Gluon propagator and comparison with the most recent lat¬ 
tice data 

From the refined Gribov-Zwanziger action, eq. (|61l) . one can immediately evaluate the tree level 
gluon propagator in linear covariant gauges, given by the following expression 


(A‘Jk)Ai(-k)) = S 


ab 


k 2 + m 2 


( k 2 + rh 2 ){k 2 + M 2 ) + 2g 2 ^ A N 


, kfjtv 

dfJ,u k 2 


+ 


Oi kjj kj, 

k 2 k 2 


( 62 ) 


A few comments are now in order. First, the longitudinal sector is not affected by the restriction 
to the Gribov region, i.e. the longitudinal component of the propagator is the same as the 
perturbative one. This is an expected result, since the Gribov region Plcg f° r linear covariant 
gauges does not impose any restriction to the longitudinal component . Second, in the limit 
a —>■ 0, the gluon propagator coincides precisely with the known result in Landau gauge |5]. In 
particular, this implies that all features of the RGZ framework derived in Landau gauge remains 
true for the transverse component of the correlation function (1621) . 

8.1 Lattice results 

As we have already mentioned before, unlike the Landau, Coulomb and maximal Abelian gauge, 
the linear covariant gauges do not exhibit a minimizing functional, a feature which is the source 
of several complications in order to construct a lattice formulation of these gauges. The study of 
the linear covariant gauges through lattice numerical simulations represents a big challenge. The 
first attempt to implement these gauges on the lattice was undertaken by [371 Ell 39], [40). More 
recently, the authors mi m sa hi sa ns] have been able to implement the linear covariant 
gauges on the lattice by means of a different procedure. 

With respect to the most recent data obtained on bigger lattices, our 

results are in very good qualitative agreement: The tree level transverse gluon propagator 
does not depend on a and, therefore, behaves like the gluon propagator in the refined Gribov- 
Zwanziger framework in the Landau gauge. On the other hand, the longitudinal form factor 
Di(k 2 ) defined by 

r) ab k k 

Oi(* 2 ) = I2 (Al(k)Al(-k)) = a . (63) 

is equal to the gauge parameter a, being not affected by the restriction to the Gribov region 
^lcg- These results are in complete agreement with the numerical data of gTl[42l[43l FI. IT). ~T0|. 
Although many properties of the Gribov region Olcg i n linear covariant gauges need to be 
further established, the qualitative agreement of our results on the gluon propagator with the 
recent numerical are certainly reassuring, providing a good support for the introduction of the 
region IIlcg- 

To end this section we point out that, recently, results for the gluon and ghost propagators in 
linear covariant gauges have been obtained through the use of the Dyson-Schwinger equations 
by 0ZlliS]. 
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9 Conclusions 


In this work we have presented a proposal for a region Ulcg in held space which is free from 
infinitesimal Gribov copies in linear covariant gauges. A local implementation of the restriction 
of the domain of integration to the region Ulcg has been worked out, resulting in the refined 
Gribov-Zwanziger action given in eqs. m, ©. 

The tree-level gluon propagator stemming from the refined action (16111 has been evaluated, being 
in good qualitative agreement with the most recent numerical lattice simulations mmmm 
Rf5] 46|. Although many geometrical properties of the region Ulcg remain to be established, we 
regard the agreement with the lattice data as a first encouraging step towards a non-perturbative 
analytic formulation of the linear covariant gauges. 

Also, we have proven that at one-loop order, the Gribov parameter 7 is independent from the 
gauge parameter a. This result is a kind of consistency check of the viability of our formulation, 
since 7 enters explicitly the computation of physical quantities like the masses of the glueballs. 

We underline that, within the class of the covariant renormalizable gauges, the non-perturbative 
treatment of Yang-Mills theories, both at analytical and numerical level, is usually carried out 
in the Landau gauge, for which a large amount of non-perturbative results has been already 
achieved. In this sense, the non-perturbative study of the linear covariant gauges is at the 
beginning. From this perspective, our results can be seen as a first step in order to address the 
Gribov issue in this class of gauges. We hope that it might be useful for the well succeeded 
interplay between numerical and analytical analysis. 

Let us end by giving a partial list of aspects of the linear covariant gauges which we intend 
to investigate in the near future. First, we should study the renormalizability properties of 
the Gribov-Zwanziger action in linear covariant gauges |26] , Second, a better understanding 
of the geometrical properties of the region Ulcg needs to be achieved. Third, the inclusion of 
quark matter fields along the lines outlined in [ 39lj can provide information about the quark 
propagator and its non-perturbative mass function in these gauges. Fourth, a detailed study 
of the Faddeev-Popov ghost propagator remains to be worked out. Fifth, the evaluation of the 
effective potential for the dimension two operators f d 4 xA T A T and f d 4 x((pip — i joj) will enable 
us to extend to higher orders the a-independence of the massive parameters ( 7 , rh, M ) entering 
the refined Gribov-Zwanziger action in the linear covariant gauges. We hope to report soon on 
those interesting topics. 
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A Conventions and notation 


In this paper, we consider pure Yang-Mills theories in four Euclidean dimensions with SU(N) 
gauge group. We choose the linear covariant gauges as our gauge condition, i.e. 

d = ab a , (64) 

where A “ denotes the gauge field, a is a gauge parameter and b a is the Lautrup-Nakanishi 
Lagrange multiplier. 

Therefore, the gauge fixed Yang-Mills action is 


Sym + S g { 


1 

4 

1 

4 


d 4 x F a F a 

Uj 1 fiu- 1 fiu 

d 4 x F a F a 

Uj 1 flu 1 LLV 


■I 


+ s I d A x <f 


J 


+ / d A x 


(m; - |i>“) 

(Mi - |i>“) + <?8„Dfc b 


(65) 


with F^ w = — d u A a ^ + gf abc A b ^A^, being the field strength, c a and c a the Faddeev-Popov 

antighost and ghost respectively, s the nilpotent BRST operator and = Y' fe <9 /t — gf abc A c ^ 
the covariant derivative in the adjoint representation. The BRST transformations for the fields 
{A, c, c, b ) are 


Su 4“ 

= -D°; b c b 

/t 


a 

_ 9 rabc b 

~ 2 J 

SC 

Sif 

II 

c s> 

sb a 

= 0. 


( 66 ) 


The gauge condition (1641) imposes a constraint on the longitudinal component of the gauge field, 
since the divergence of the transverse sector is identically vanishing. Therefore, it is useful to 
decompose the gauge field into transverse and longitudinal components, 


A Ta — (5 — JS& 

A* — l °l‘- v Q2 J 

a La _ . a 

/t Q2 v ' 


(67) 


Clearly, the choice a = 0 in ea. ffMl) makes the gauge condition equivalent to the Landau gauge, 
which implies the vanishing of the longitudinal component of the gauge field. 
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